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1. INTRODUCTION 

Electrical transformers the facilities in electrical power networks that the correct operation of it has 
played main role in providing reliability of these systems. By noticing the costs of transformers and 
permanent connection of transformers to transmission and distribution systems, any action to increase the life 
cause improvement in the operation of power systems. Life of the electrical transformer and loading 
capabilities of this facilities depend on several parameters. Oil and paper in transformers are two insulation 
materials, that the integrity of them warranties the life of transformers. Hot spot temperature (HST) and top 
oil temperature (TOT) are two criteria to evaluate the integrity of these two insulation materials and 
transformers as important facility in power network [1]-[3]. These criteria determine the temperature 
limitation of transformer for loading capability, for example the IEFE standard suggests that the TOT rise 
over ambient temperature 65°C and depending on it has considered 110°C for HST [4]-[6]. It means that 
increasing the temperature of transformer over this standard value causes the degradation of insulation life of 
transformer. According to research, from 90°C to 110°C, the aging rate of tensile strength approximately 
doubles every 5°C to 10°C temperature rise [2], [3], [7]-[9]. A. A. Taheri, A. Abdali and A. Rabiee [10] a 
thermal resistance with solar radiation is used to explain the heat transfer theory in transformer. The effect of 
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long term ambient temperature was investigated in [11], [12]. There are some dynamic factors in oil 
immersed transformer such as ambient temperatures, winding losses and oil viscosity of transformer [3], 
[13]. In this article, two thermal models, classic thermal model as famous and reliable thermal model and 
dynamic thermal model as improved thermal model, has been used. Moreover, TOT and HST is calculated 
for both thermal model by IEEE standard method that has been available in [3] and also these temperatures 
are calculated by proposed method. To validate and show the tolerance of numerical technique, the results of 
models compared with measured value. 


2. LITERATURE REVIEW 
2.1. Thermal model 
2.1.1. Classic thermal model 

Thermal model is used to estimate the critical temperature of the transformer. The classic one was 
given base on the heat generated and distributed in the windings and oil of transformer [3], [9]. Figure 1 
shows temperature distribution in the transformer. 
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Figure 1. The temperature distribution in the transformer 


The classic thermal model is used by a differential equation as [3]: 


dôtop 








To i —Frop + Ou + Oamv (1) 
The solution as: 
= -© 
Prop = (Ou F Oamv 7 topi) (1 =e o ) + topi (2) 
and 
I2uR+1]” 
Oy = Op | a | (3) 
co 
T=% (4) 
fl 


where: 8;9p=Top-oil rise over ambient temperature (C) 
6,,=Ultimate top-oil rise for load L (C) 
6;=Initial top-oil rise for t=0 (C) 
6,,=Top-oil rise over ambient temperature at rated load 
To=Time constant 
n=Oil exponent 
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I=Specified load 

L-atea=Rated load 

Ipu=Ratio of specified load to rated load 
The transient rise of the hottest spot of the winding higher than the upper limit of the oil temperature is 
determined by the following formula [3]: 


=E 


40y = (40y u — 40y) £ — expTo | +4ðyi (5) 


That the initial value for hot-spot temperature rise over top-oil temperature is given by: 


AOy i = AOy pl? (6) 


ipu 
The maximum hot spots are as follows: 

Abn uy = AO x, rl pu (7) 
The nominal hottest temperature spot value is given by: 

AOR = AOu AR = AGOro,R (8) 


where: A@;,=Hottest temperature spot rise over top oil temperature (C) 

A@y y=Maximum hottest temperature spot rise over maximum oil temperature for load(C) 

A@,, ;=Initial hottest temperature spot rise over top oil temperature for t = 0 (C) 

To=time constant of the transformer winding (h) 

t=Loading time (h) 

40y p=Nominal hot spot rising above maximum oil temperature (C) 

I=Current of load per unit 

M=The exponent determines the non-linearity as shown in Table 1 

The manufacturer has specified the amount of n for each mode of cooling. The amounts of n and m 
for each cooling method are shown in Table 1 [3]. 


Table 1. Exponents for temperature determination formulas 











cooling M N 
OA 0.8 0.8 
FA 0.8 0.9 
Non-directed FOA or FOW 0.8 0.9 
Directed FOA or FOW 1 1 
HST will be calculated by (9). 
Gust = AOror + APusr + Oamb (9) 


2.2.2. Dynamic thermal model 

Firstly, A dynamic thermal model of electrical transformer has been offered by Swift based on the 
principals of heat transfer theory [14]-[17] after that an improved model has been introduced by Susa by 
assuming the non-linear thermal oil resistance base on swift’s approach. In Susa’s model the variations of 
viscosity with temperature were considered [9]. TOT of dynamic thermal model is governed by (1). 


1+RXK? dOoit 


6,:)-8. 1+n 
sae X Mou X Ofi = Hpu X To X eS (10) 


dt OF, 





Oil viscosity is added as a parameter. This can be determined by [18]. 
í 2797.3 )-( 2797.3 | 
Upu =e Otop+t273 8 f1+273 ad 1) 


where: K=load current per unit 
R=proportionality of load to no-load losses 
Tg=top oil time constant, (min) 
0,,=rated top oil rise over ambient, (K) 
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O i=top oil temperature gradient, (K) 

6 amp =ambient temperature, (K) 

Hpy=oil viscosity per unit 

n=empirical constant, the amount of n has been specified for each cooling mode in Table 2 [2]. 


Tabel 2. Experimental value of n 








Oil circulation n 
With cooling without cooling 
Initial oil circulation, Speed=0, (cold start) ONAF/OFAF ONAN ONAN 
0.5 0 0.25 
Initial oil circulation, Speed>0, (transformer on load) 0.2 0.25 0.25 





HST of dynamic thermal model is given as [19]: 





2 = dðns , ((Ons—9oit)"** 
[K x Pediat x Upu x AGns rated = Upu x Twd,rated x dt F ( Aor iad (12) 
s,rate 


Pwd,pu, is the load loss’s dependence on temperature that can be determined by (13) 


9nst Ox 9nsratedt Ox 
Pwapu = Pacpu X (SS) +P, x (tsrereattr 13 
wd,pu dc,pu OnsratedtOk eddy,pu OnstOx ( ) 


The behavior of the eddy losses and DC losses are shown in Peddy and Pdc respectively. DC loss is 
directly proportional to temperature, while eddy current loss is inversely proportional to temperature. Ok is 
the temperature coefficient of loss correction. Aluminum is equal to 225 and copper is equal to 235.Oil time 
constant (to) is one of the most important parameters that should be determined accurately. Note that, there 
are specific methods for each thermal model to determine to. The oil time constant for dynamic thermal 
model is given by [20]: Note that, (14) is based on minute and (15) is per hours and depends on whether per 
hours or minutes have to be changed and for data which collected per hours it shouldn’t multiply by 60, it 
means that the oil time constant does not need to convert to minutes. 


Cen—oil X49 ot 
To a th—oil oilrated x 60 14) 
P 
total 
To = Cth-oilX 49 oil rated (15) 
P 
total 


3. RESEARCH METHOD 
3.1. Solving techniques 
3.1.1. IEEE standard method 
IEEE standards C57.91-1995 has offered an exponential method based on the (14) to calculate the 
TOT to determine the HST of transformers for classic thermal model [3]. In (14) is derived from (1) [21]. 


-1 
48ro = (A8roy — 48ro;) |1 — exp] + Ab ro, (16) 


The ultimate value of top oil acsent for each stage is calculated by, 


n 
CpuuR+1) 
(R+1) 


Grou = fi | (17) 
For calculating the initial value for the first time, t=0, knowing the record of loading profile can be 

so useful. The current at t=0, determined by RMS value of the loading currents for the former loadings [3]. 

Equation (3.5) can be used to estimate the initial value for other steps to calculate the top-oil temperature rise. 


a 18) 


— g, [Cui 
48ro,i = bri | R 


3.1.2. Numerical method 


Numerical methods are used to solve the ordinary differential equation. In this study, two types of 
numerical solvers are compared to determine the thermal-electrical models. The first model is related to 
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Euler's method and is the simplest solver. Although this model is not an appropriate solver in most 
utilization, it reveals some common characteristics among different types of solvers. Runge-Kutta is the 
second type of solver, which is used in this paper. The fourth-order RungeKutta technique is the most widely 
used solution method. This solution technique, which begins with an initial value (to,yo), has four slopes, 


h h 
K, = f Cto; Yo), K2 = f (to +5 Yo +3 K1) 


= (19) 
K, =f (ty +i, y +K), Ka = f (to + h, Yo + hK3) 
With these slopes, the next value of the dependent variable can be determined as, 
K1+2K2+2K3 +K. 
Yk = Yk-1 + h= (20) 


TOT and HST of transformer can be computed for dynamic thermal model according to the algorithm 
has been given in Figure 2. The input data can be classified in two groups, constant and variable input data. 
Constant data are n, m, weight of the oil per kg, load losses and no-load losses, oil time constant and winding 
time constant, DC and eddy losses, And the hot spot temperature under rated load. While variable input data 
that can be changed by time include current or loading of transformers and ambient temperature. 


Input n, m, Opm 


» -m_oil. Oa Py.Py.Tust. O, 


topoil _ initial >° 


K=load current, Pac-pu-PED-pu- OgrsT _iniriar > Pusr 


_rate 


PaA 
Calculate R= =— „and T, by (14. 15) 
N_L 








Calculate 4p, by 1) 


Calculate TOT, equation (10) , by 


Euler and RK4 method, algorithm 
@) and 6) 


YES : 








Calculate M, by aD 


pu 


Calculate P, for t=0 by (13) 


wd pu 


Calculate HST, equation (12) , by 


Euler and RK4 method, algorithm 
®© and G) 








> 


Plot TOT, HST 


Figure 2. HST and TOT for dynamic thermal model 
*: N is total number of load current that are collected 
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3.2. The loss of life of the transformer 

The quality of paper and oil which are applied in electrical transformer is directly proportional to the 
efficiency of it. However, the deterioration of paper can determine the loss of life (LOL) the mineral oil 
immersed transformers. Heat, oxygen and moisture are the most important factors to deteriorate the paper of 
electrical transformers. Insulation degradation can be simulated as a quantity per unit and calculated as [4], 
[22]; 


B 
per — unit — life = Aclrnsr (21) 


where: A=is an adjusted constant based on the temperature determined for a per unit life, 
B=is the slope rate of aging. 
IEEE standard C57.91-1995 considered B=15000 is appropriate and is used in per-unit-life equation 
for loading of transformers. For @457=110 C, and per-unit-life=1, the value of A constant becomes 
9.8x10-18. The aging factor can be calculated [3], [23] 








15000 15000 l 
Fy =e 383 0H+273 (22) 

The equivalent aging of the electrical transformer based on reference temperature for hottest spot 
is 110°Ccould be calculated as [3], [24], [25]: 


N 
Fega = iy (23) 
where: Feoa=aging factor 
n=time interval index 
N=total time duration 
At,=time interval 
And, the percentage of LOL can be determined as [26], [27]: 


Frgaxtx100 
NIL 


%loss — of — life = (24) 
3.3. Emprical tests for case of study 
3.3.1. Three phase 2500 KVA transformer 

The ratio of voltage of 2500 KVA transformer is 20/0.4 KV, and without external cooling, ONAN. 
Three different loading currents are considered in this paper that are shown in Table 3. For this transformer 
the low voltage and high voltage windings have one and two axial cooling duct respectively. The transformer 
equipped by 17 thermocouples located as follows: One is mounted 50 mm under the tank cover of electrical 
transformer, 2 sensors in the oil enclosure,4 sensors is located on the bottom of outside surface of the tank, 2 
in the bottom of the tank, 4 in duct outlet in top and bottom, 4 in duct inlet in top and bottom. 


Table 3. Load levels for transformer 








Time intervals (minutes) Load (PU) 
0-300 1 
300-382 2 
382-415 0 





3.3.2. Three phase 30 MVA transformer 

The voltage ratio of this transformer is 33/11 KV. Also, cooling system is oil natural and air force 
(ONAP). 3 thermometers are used to indicate the winding and top oil temperatures for high and low voltage 
windings. Also, they are used for collecting ambient temperature. This data is necessary to verify the 
suggested techniques provided in this article. The output of the thermometer is a function of the 4 to 20 mA 
current that must be converted to temperature (C). 


4. RESULTS AND DISCUSSION 
HST of each transformer has been plotted for both distribution transformers. LOL of transformers 
for each thermal model has been compared with measured value and were indicated in Tables 4 and 7 for 
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classic thermal model and Tables 5 and 8 for dynamic thermal model. Finally, Tables 6, and 9 tabulate the 
error of numerical methods rather than exponential method. By comparing the error of each thermal model 
for 2500 KVA transformer, which is demonstrated in Tables 4 and 5, the less error of dynamic thermal model 
rather than classic thermal model is comprehensible. These errors for numerical methods are less than 
exponential method. 


4.1. 2500 KVA transformer 

2500 KVA transformer was tested for 415 minutes. Figure 3 shows the comparison between the 
measured HST and HST achieved by the numerical methods and exponential method. Graphs have been 
plotted based on the classic model. Figure 4 shows the comparison of the hot spot temperature of a 2500 
KVA transformer with the hot spot temperature calculated using two numerical methods (RungeKutta and 
Euler methods).This graph has been plotted based on the dynamic thermal model. 
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Figure 3. Hot spot temperature of 2500 KVA transformer: Classic thermal model 


HST for 2500KVA transformer 
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Figure 4. Hot spot temperature of 2500 KVA transformer: Dynamic thermal model 


In Tables 4 and 5 show the comparison between the errors achieved for classis and dynamic thermal 
models respectively. These errors given the differences between the numerical methods and exponential 
methods with measured. In addition, the loss of life of the 2500 KVA transformer is compared with measured 
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value in the Table 6 for two types of the thermal models. Thermal models have been calculated based on R-K 
technique. This comparison shows that the accuracy of the dynamic thermal model is higher than that of the 
classical thermal model. 


Table 4. Classic thermal model of 2500 KVA transformer (ERROR OF LOL) 
Error for RK4 |(LOLNUM — LOLMeasured)| Error of exponential|(LOLIEEE — LOLMeasured)| 
1.7052 2.7643 











Table 5. Dynamic thermal model of 2500 KVA transformer (ERROR OF LOL) 
Error for RK4 method |(lol[RK4 — lolMeasured)| Error for Euler method |(lolEuler — lolMeasured)| 
0.7806 1.4912 











Table 6. Error of thermal models: Different solving technique 








numerical technique numerical technique exponential technique 
(Dynamic thermal model ) (Classic thermal model) (Classic thermal model) 
0.2358 2.2503 2.7643 





5. CONCLUSION 

Several thermal models and solution are suggested for solving various thermal models to calculate 
the TOT and HST of transformer due to loss of life evaluation. These methods involved are linear and 
nonlinear regression and exponential methods that are proposed by IEEE standard. Numerical methods can 
be used as solution to calculate the thermal model, R-K technique is an improved numerical method proposed 
in this research for solving two types of thermal models the first one is classical thermal model that was 
suggested by IEEE standard C57.91-1995 and the second one is thermal-electrical model or dynamic thermal 
model that assumed number of variable parameters such as viscosity of oil, ambient temperature, eddy losses. 
According to the results the proposed technique for solving the thermal models shows improved accuracy 
over the traditional methods. 
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